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AinendedClaim4 to Obtain Allowance Qf Case 

4. (Cuirently Amended) A method for the generation of a three dimensional fiactal 
" subsurface structure by Voronoi tessellation and computation of gravity response of such fractal 
structure, by generation of fractal subsurface structures and computation of a forward gravity 
response of such structure for delineation and modeling of an underlying anomalous object, said 
method comprisixjg: 

(a) selecting Voronoi centers at a pliurality of locations over a region of interest, the 
Voronoi centers being represented by x, y, z co-ordinates; 

(b) generating an initial model of a subsurface fractal geological object, having variation 
'i m^a physical property in lateral and vertical directions; the initial nipdel being generated by 

tessellating the Voronoi centers, and assigning values of physical property variations during 
generation of the model on the basis of pre-determined assumptions in the model; 

(c) assigning different discernable representations to regioM 
'^^ djbfferjsnt physical propertie to enable demarcation; 

wherein the fractal subsurface model is generated bv a modified Voronoi tessellation 
technique which comprises modifying the Voronoi tessellation bv using nonn. where LF= (x- 
qiV^^ , where x is an arbitrary point, qj^is point whose distance has to be calculated, and o is an 
exponent which can assume finctional real numbe rs and negative real numbers therebv enabling 
greater range ftur generation of different fractal sub-surfaces, which are closer to a real gdological 
situation, 

wherein the computation of a gravity response due to a fractal subsurface generated by 
[ modified Voronoi tesselation is carried out by an analytical expression which comprises: 
(a) demarcating boundaries of regions having a different physical property in the 



PAGE 2/4'RCVDAT8/14/200712:13:48PM [Eastern Daylight 



NIXON S VANDERHYEPC Fax:703-816-4100 



Aus 14 2007 12:20 



P. 03 



I tessellated region, the boundaries forming a polygonal shape in 2-dimensional space; 
V = Gp ±{W arccos{{x, lr^){x,^, lr,,,)+{y, /r,)fe,, Ir,,,)} 



- arc sin — ^ ^ + arc sin 



Wdere 8== +1 if pi is positive, S= -1 if pi is negative, 

W= +1 if nii is positive, W= -1 if mj is negative, 

*Z' is depth aad 'n' is number of sides in the polygon. 
G is universal gravitational constant, p is the density of the tessellated regions; 

yi~ yM ^ ^i^^M 

^ ^ ^i-^M . yj-y^ yj 

Hi ~ \ ^ > .... 

^UM ^hM .... 

. , 

Pi = perpendicular to a side of an irregular geometrical body from a point at which 
i_ aQpinaly is being calculated; 

qi ^ <^os(6j)> where 0^ is an angle between a side of the irregular body and line 
. Joining the first point (A) of the arm from the point at which anomaly is being computed; 



,1/2-^ 
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fi = cos((pi), where 91 is an angle between a side of the irregular body and line 
" joining the end point (B) of the arm from the point at which anomaly is being computed, 
mi ^ angle subtended by an arm of the irregular body at the point where anomaly 
is being computed; 

ri = distance of a point at the boundary of irregular body from the point at which 
. anomaly is being computed; 

ri+i = distance of next consecutive point at the boundary of the irregular body from 
the point at which anomaly is being calculated; 

^U+i ^ length of the side of the irregular body; 
. i = is a subscript, which is used to select the next consecutive point on the surface 

of the body; 

where the effect of a common arm of the adjacent polygon is removed; 
" (b) repeating the gravity response computation due to the polygon for all adjacexxt 

polygons of diffident physical properties using the demarcated polygon boundary; 

(c) adopting the process of gravity response computation for tesseUated regions lying at 
different depths to obtain a response; 

(d) integrating the response using Simpson/Gauss quadratvire formula at plurality of grid 
nodes overlain on the region of interest; and 

(e) creatingamodelof the underlying anomalous object. 
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